Gaussian measures of entanglement versus negativities: 
the ordering of two-mode Gaussian states 



Gerardo Adesso and Fabrizio Illuminati 

Dipartimento di Fisica "E, R. Caianiello ", University di Salerno; CNR-Coherentia, Gruppo di Salerno; 
and IN FN Sezione di Napoli-Gruppo Collegato di Salerno, Via S. Allende, 84081 Baronissi (SA), Italy 

(Dated: June 20, 2005) 

In this work we study the entanglement of general (pure or mixed) two-mode Gaussian states of continuous 
variable systems by comparing the two available classes of computable measures of entanglement: entropy- 
inspired Gaussian convex-roof measures, and PPT-inspired measures (negativity and logarithmic negativity). 
We first review the formalism of Gaussian measures of entanglement, adopting the framework introduced in 
[M. M. Wolf et al, Phys. Rev. A 69, 052320 (2004)], where the Gaussian entanglement of formation was de- 
fined. We compute explicitely Gaussian measures of entanglement for two important families of nonsymmetric 
two-mode Gaussian states, namely the states of extremal (maximal and minimal) negativities at fixed global 
and local purities, introduced in [G. Adesso et al., Phys. Rev. Lett. 92, 087901 (2004)]. This analysis allows to 
compare the different orderings induced on the set of entangled two-mode Gaussian states by the negativities 
and by the Gaussian measures of entanglement. We find that in a certain range of values of the global and local 
purities (characterizing the covariance matrix of the corresponding extremal states), states of minimum neg- 
ativity can have more Gaussian entanglement of formation than states of maximum negativity. Consequently, 
Gaussian measures and negativities are definitely inequivalent measures of entanglement on nonsymmetric two- 
mode Gaussian states, even when restricted to a class of extremal states. On the other hand, the two families of 
entanglement measures are completely equivalent on symmetric states, for which the Gaussian entanglement of 
formation coincides with the true entanglement of formation. Finally, we show that the inequivalence between 
the two families of continuous-variable entanglement measures is somehow limited. Namely, we rigorously 
prove that, at fixed negativities, the Gaussian measures of entanglement are bounded from below. Moreover, we 
provide some strong evidence suggesting that they are as well bounded from above. 

PACS numbers: 03.67.Mn, 03.65.Ud 



I. INTRODUCTION 

Quantum information with continuous variables (CV) Q0I 
is a flourishing field dedicated to the manipulation of the in- 
formation using quantum states governed by the laws of quan- 
tum mechanics. This approach contrasts with the usual meth- 
ods involving discrete-spectrum observables (such as, e. g., 
polarization, spin, energy level) of single photons, atoms or 
ions. The ability of quantum states with continuous spectra to 
implement quantum cryptography |3|], quantum teleportation 
Hp |a Q entanglement swapping [9|], dense coding 
llOtl . quantum state storage flUl . and, to some extent, quan- 
tum computation lll2ll processes, brings up new and exciting 
perspectives. 

The crucial resource enabling a better-than-classical manip- 
ulation and processing of information is CV entanglement, 
introduced for the first time in the landmark paper by Ein- 
stein, Podolski and Rosen IQI in 1935. There, it was shown 
that the simultaneous eigenstate of relative position and to- 
tal momentum of two particles (or of a two modes of the ra- 
diation field) contains perfect quantum correlations, i.e. infi- 
nite CV entanglement. While this state is clearly an unphys- 
ical, unnormalizable state, it can be approximated arbitrarily 
well by two-mode squeezed Gaussian states with large enough 
squeezing parameter. The special class of Gaussian states 
(which includes thermal, coherent, and squeezed states), thus 
emerges quite naturally in the CV scenario. These entangled 
states can be easily produced and manipulated experimen- 
tally, and moreover their mathematical description is greatly 
simplified due to the fact that, while still living in a infinite- 



dimensional Hilbert space, their relevant properties (such as 
entanglement and mixedness) are completely determined by 
the finite-dimensional covariance matrix of two-point corre- 
lations between the canonically conjugated quadrature opera- 
tors. Therefore, clarifying the characterization and quantifica- 
tion of CV entanglement in two-mode and, eventually, mul- 
timode Gaussian states stands as a major issue in the field of 
C V quantum information, as the amount of entanglement con- 
tained in a certain state directly quantifies its usefulness for 
information and communication tasks like teleportation |8|. 

For the prototypical entangled states of a CV system, the 
two-mode Gaussian states, much is known about entangle- 
ment qualification, as the separability is completely character- 
ized by the necessary and sufficient PPT criterion (positivity 
of the partially transposed state) i\M \ and also with regard 
to its quantifcation. Concerning the latter aspect, the negativ- 
ity (quantifying the violation of the necessary and sufficient 
PPT condition for separability) is computable for all two- 
mode Gaussian states Ml 511 . Moreover, for symmetric two- 
mode Gaussian states also the entanglement of formation is 
computable [16], and it turns out to be completely equivalent 
to the negativity for these states. 

Another measure of CV entanglement, adapted for the class 
of Gaussian states, has been introduced in Ref. 11711 . where 
the Gaussian entanglement of formation (an upper bound to 
the true entanglement of formation) was defined as the cost 
of producing an entangled mixed state out of an ensemble of 
pure, Gaussian states. While the Gaussian entanglement of 
formation coincides with the true one for symmetric states, at 
present it is not known whether this equality holds for non- 
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symmetric states as well El. 

In this work, aimed at sheding new light on the quantifica- 
tion of entanglement in two-mode Gaussian states, we com- 
pute the Gaussian entanglement of formation and, in general, 
the family of Gaussian entanglement measures, for two differ- 
ent classes of two-mode Gaussian states, namely the states of 
extremal, maximal and minimal, negativities at fixed global 
and local purities |19, 20]. We find that the two families of 
entanglement measures (negativities and Gaussian measures) 
are not equivalent for nonsymmetric states. Remarkably, they 
may induce a completely different ordering on the set of en- 
tangled two-mode Gaussian state: a nonsymmetric state qa 
can be more entangled than another state qb, with respect to 
negativities, and less entangled than the same state qb, with 
respect to Gaussian measures of entanglement. However, the 
inequivalence between the two families of measures is some- 
how bounded: we show that, at fixed negativities, the Gaus- 
sian entanglement measures are rigorously bounded from be- 
low. Moreover, we provide strong evidence hinting that they 
should be bounded from above as well. 

The paper is organized as follows. In Sec. |II]we set up 
the notation and review the basic properties of Gaussian states 
of CV systems. In Sec. [|I1]we review the main results on the 
characterization of separability in Gaussian states, introducing 
also two families of measures of entanglement, respectively 
the negativities and the Gaussian entanglement measures. In 
Sec. HVI we compute the latter for two-mode Gaussian states, 
solving the problem explicitely for the states of extremal nega- 
tivities at fixed purities, described in Sec. IIV Al In Sec. fVlwe 
compare the orderings induced by negativities and Gaussian 
measures on the set of extremal two-mode Gaussian states. In 
Sec. IVII we compare the two families of measures for generic 
two-mode Gaussian states, finding lower and upper bounds 
on one of them, when keeping the other fixed. Finally, in Sec. 
IVIII we summarize our results and discuss future perspectives. 

H. GAUSSIAN STATES: DEFINITIONS AND NOTATION 

A continuous variable (CV) system is described by a 
Hilbert space H — <S)2=i Hi resulting from the tensor product 
structure of infinite dimensional Fock spaces Hi's. Let at be 
the annihilation operator acting on Hi, and q\ = (a, + aj) and 
f>i = (a,i — a\) I 'i be the related quadrature phase operators. 
The corresponding phase space variables will be denoted by 
qi and p^. Let X = (qi,pi, ■ ■ ■ , q n ,Pn) denote the vector of 
the operators fa and pi. The canonical commutation relations 
for the Xi can be expressed in terms of the symplectic form 57 

[Xi, Xj] = 2iflij , 



with £1 = u> , lu = 

i=l 




The states of a CV system can be equivalently described 
by a positive trace-class operator (the density matrix q) or 
by quasi-probability distributions such as the Wigner function 



ll2lll . States with Gaussian characteristic functions and quasi- 
probability distributions are referred to as Gaussian states. 
Such states are at the heart of information processing in CV 
systems [2] and are the subject of our analysis. By definition, 
a Gaussian state q is completely characterized by the first and 
second statistical moments of the quadrature field operators, 
which will be denoted, respectively, by the vector of first mo- 
ments X = (iXx), (X\), . . . , (X n ), (X n )^j and the covari- 
ance matrix (CM) er of elements 

an = ~ (X t X 3 + X 3 Xi) - (Xi) (Xj) , (1) 

where, for any observable 6, the expectation value (6) = 
Tr (go). Notice that the entries of the CM can be expressed as 
energies by multiplying them by the quantity hw, where w is 
the frequency of the considered mode. In fact, for any ?i-mode 
state (even non Gaussian) the quantity huTi er /4 is simply the 
average of the non interacting Hamiltonian Yl7=i ( a ! a i+l/2). 
First moments can be arbitrarily adjusted by local unitary op- 
erations (displacements), which cannot affect any property 
related to entropy or entanglement. Therefore, they will be 
unimportant to the present scope and we will set them to in 
the following, without any loss of generality. 

The canonical commutation relations and the positivity of 
the density matrix g imply 

er + id > , (2) 

Inequality (|2ji is the necessary and sufficient constraint the ma- 
trix a has to fulfill to be a CM corresponding to a physical 
Gaussian state 1221 12311 . More in general, the previous condi- 
tion is necessary for the CM of any, generally non Gaussian, 
state. We note that such a constraint implies er > 0. 

A major role in the theoretical and experimental manipula- 
tion of Gaussian states is played by unitary operations which 
preserve the Gaussian character of the states on which they 
act. Such operations are all those generated by Hamiltonian 
terms at most quadratic in the field operators. As a conse- 
quence of the Stone- Von Neumann theorem, any such uni- 
tary operation at the Hilbert space level corresponds, in phase 
space, to a symplectic transformation, i.e. to a linear trans- 
formation S which preserves the symplectic form f2, so that 
n = S T nS. Symplectic transformations on a 2n-dimensional 
phase space form the (real) symplectic group Sp(2n,R)- Such 
transformations act linearly on first moments and by con- 
gruences on covariance matrices: er i— * S T crS. One has 
DetS 1 = 1, V5 € Sp(2n,B.)- Ideal beam splitters, phase 
shifters and squeezers are all described by some kind of sym- 
plectic transformation. A particularly important symplectic 
transformation is the one realizing the decomposition of a 
Gaussian state in normal modes. Th roug h this decomposi- 
tion, thanks to Williamson theorem 12411 . the CM of a n- 
mode Gaussian state can always be written in the so-called 
Williamson normal, or diagonal form 

cr = S T vS , (3) 
where S € <Sp(2n,B) an d v i s me CM 

v = diag(z/i , . . . , v n , v n ) , (4) 
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corresponding to a tensor product of thermal states with a di- 
agonal density matrix g® given by 



(5) 



where denotes the number state of order k in the Fock 
space Hi . 

The quantities fi's form the symplectic spectrum of the CM 
er, and they can be computed as the eigenvalues of the matrix 
|ifier|. Such eigenvalues are in fact invariant under the action 
of symplectic transformations on the matrix er. The symplec- 
tic eigenvalues vi encode essential informations on the Gaus- 
sian state er and provide powerful, simple ways to express its 
fundamental properties. For instance, in terms of the symplec- 
tic eigenvalues Ui, the uncertainty relation reads 



Vi > 1 . 



(6) 



Moreover, the entropic quantities of Gaussian states can be 
as well expressed in terms of their symplectic eigenvalues and 
invariants 1 20] . Notably, the purity Tr g 2 of a Gaussian state g 



is simply given by the symplectic invariant Dct er 
being |25] 



n; 



— Tr g — 



1 



VDct< 



(7) 



A. Two-mode states 

This work is focused on two-mode Gaussian states: we 
thus briefly review here some of their basic properties. The 
expression of the two-mode CM er in terms of the three 2x2 
matrices ct, f3, 7, that will be useful in the following, takes 
the form 



a. 7 
7 T 



(8) 



For any two-mode CM er there is a local symplectic operation 
Si = 5j © 5*2 which brings er in the so called standard form 



Sftr Si — er s ^ 



a c+ 

a c_ 

c+ b 

V c_ b 



(9) 



States whose standard form fulfills a = b are said to be sym- 
metric. Let us recall that any pure state (/i = 1) is symmetric 
and fulfills c + = — c_ = y/a 2 — 1. The correlations a, b, 
c+, and c_ are determined by the four local symplectic invari- 
ants Deter = (ab — c 2 + ){ab — c 2 .), Detct = a 2 , Det/3 = b 2 , 
Det7 = c + c_. Therefore, the standard form corresponding 
to any CM is unique (up to a common sign flip in c_ and c + ). 

For two-mode states, the uncertainty principle Ineq. (0 can 
be recast as a constraint on the Spyj) invariants Deter and 
A(cr) = Deta + Det/3 + 2 Det7 Q: 



A(er) < 1 + Deter 



(10) 



The symplectic eigenvalues of a two-mode Gaussian state 
will be denoted as v_ and v + , with v_ < v + , with the uncer- 
tainty relation (|6} reducing to 



v- > 1 



(ID 



A simple expression for the v T can be found in terms of 
the two Sp( 4 r) invariants (invariants under global, two-mode 

symplectic operations) din 



A(er) =F V A2 <» -4 Deter. 



(12) 



III. ENTANGLEMENT OF GAUSSIAN STATES 

In this section we recall the main results on the qualification 
and quantification of entanglement for Gaussian states of CV 
systems. 

A. Qualification: PPT criterion 

The positivity of the partially transposed state (Peres- 
Horodecki PPT criterion |28|) is necessary and sufficient for 
the separability of two-mode Gaussian states 11411 and, more 
generally, of all (1 + n)-mode Gaussian states under 1 x n- 
mode bipartitions 12911 and of symmetric and bisymmetric 
(m + nV-mode Gaussian states under m x n-mode biparti- 
tions 13011 . In general, the partial transposition g of a bipar- 
tite quantum state g is defined as the result of the transposi- 
tion performed on only one of the two subsystems in some 
given basis. In phase space, the action of partial transposi- 
tion amounts to a mirror reflection of one of the four canoni- 
cal variables 11411 . The CM er is then transformed into a new 
matrix er which differs from er by a sign flip in Det 7. There- 
fore the invariant A(er) is changed into A(er) = A(er) = 
Det a. + Det (3 — 2 Det 7. Now, the symplectic eigenvalues 
f>zr of er read 



\ 



A(er) =p A/A 2 (er) -4 Deter 



(13) 



The PPT criterion for separability thus reduces to a simple 
inequality that must be satisfied by the smallest symplectic 
eigenvalue of the partially transposed state 



v- >1, 

which is equivalent to 

A(er) < Deter + 1 



(14) 



(15) 



Moreover, the above inequalities imply Det 7 = c + c_ < 
as a necessary condition for a two-mode Gaussian state to be 
entangled. Therefore, the quantity i>_ encodes all the qualita- 
tive characterization of the entanglement for arbitrary (pure or 
mixed) two-mode Gaussian states. 
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B. Negativities 

From a quantitative point of view, a measure of entangle- 
ment which can be computed for general Gaussian states is 
provided by the negativity J\f, first introduced in Ref. [3111 _ 
later thoroughly discussed and extended in Refs. flU l32ll to 
CV systems. The negativity of a quantum state q is defined as 



(16) 



where q is the partially transposed density matrix and ||d|| i = 
Tr|6| stands for the trace norm of the hermitian operator 6. 
The quantity M(q) is equal to | A» |, the modulus of the 
sum of the negative eigenvalues of q, quantifying the extent 
to which q fails to be positive. Strictly related to M is the 
logarithmic negativity Ej^, defined as Ej\f = log ||f?||i, which 
constitutes an upper bound to the distillable entanglement of 
the quantum state q and is related to the entanglement cost 
under PPT preserving operations 13 311 . Both the negativity and 
the logarithmic negativity have been proven to be monotone 
under LOCC (local operations and classical communications) 
a crucial property for a bona fide measure of 
entanglement. Moreover, the logarithmic negativity possesses 
the nice property of being additive. 

For any two-mode Gaussian state q it is easy to show that 
both the negativity and the logarithmic negativity are simple 
decreasing functions of V- 115112011 



\q\\ = r— => N(q) = max 



0.- 



2£_ 



Em{q) = max [0, - log£>_] 



(17) 



(18) 



These expressions directly quantify the amount by which the 
necessary and sufficient PPT condition ( 1141 for separability 
is violated. The symplectic eigenvalue i>_ thus completely 
qualifies and quantifies (in terms of negativities) the entangle- 
ment of a two-mode Gaussian state er: for i>_ > 1 the state 
is separable, otherwise it is entangled. Finally, in the limit of 
vanishing i>_, the negativities grow unboundedly. 



C. Entanglement of Formation 

In the special instance of symmetric two-mode Gaussian 
states, the entanglement of formation (EoF) |35], can be com- 
puted as well II 1611 . We recall that the EoF Ep of a quantum 
state q is defined as 

E f {q)= min y>£(|^)), (19) 

where the minimum is taken over all the pure states realiza- 
tions of q: 



The asymptotic regularization of the entanglement of forma- 
tion coincides with the entanglement cost Ec(q), defined as 
the minimum number of singlets (maximally entangled anti- 
symmetric two-qubit states) which is needed to prepare the 
state q through LOCC 10. 

The optimal convex decomposition of Eq. ( 1191 has been 
found for symmetric two-mode Gaussian states, and turns 
out to be Gaussian, that is, the absolute minimum is realized 
within the set of pure two-mode Gaussian states 11611 . yielding 



Ep = max [0, h(i/-) 



(20) 



with 



h(x) 



(1 +x) 2 
4x 



log 



(1 



4,j; 



(1 - xf 

4a; 



log 



(1 



4.T 



(21) 

Such a quantity is, again, a monotonically decreasing function 
of thus providing a quantification of the entanglement of 
symmetric states equivalent to the one provided by the nega- 
tivities. 

As a consequence of this equivalence, it is tempting to con- 
jecture that there exists a unique quantification of entangle- 
ment for two-mode Gaussian states, embodied by the smallest 
symplectic eigenvalue £>_ of the partially transposed CM, and 
that the different measures simply provide trivial rescalings 
of the same unique quantification. In particular, the ordering 
induced on the set of entangled Gaussian state is uniquely de- 
fined for the subset of symmetric two-mode states, and it is 
independent of the chosen measure of entanglement. How- 
ever, regrettably, in Sec. [V]we will indeed show that different 
measures of entanglement induce, in general, different order- 
ings on the set of nonsymmetric two-mode Gaussian states. 



D. Gaussian convex-roof extended measures 

In this subsection we consider a family of entanglement 
measures exclusively defined for Gaussian states of CV sys- 
tems. The formalism of Gaussian entanglement measures 
(Gaussian EMs) has been introduced in Ref. 1 17] where the 
Gaussian EoF has been defined and analyzed. Furthermore, 
the framework developed in Ref. 11711 is general and enables 
to define generic Gaussian EMs of bipartite entanglement by 
applying the Gaussian convex roof, that is, the convex roof 
over pure Gaussian decompositions only, to any bona fide 
measure of bipartite entanglement defined for pure Gaussian 
states. The original motivation for the introduction of Gaus- 
sian EMs stems from the unfortunate fact that the optimiza- 
tion problem Eq. \\9\ for the computation of the EoF of non- 
symmetric two-mode Gaussian states has not yet been solved, 
and it stands as an open problem in the theory of entangle- 
ment 1 18]. However, the task can be somehow simplified by 
restricting to decompositions into pure Gaussian states only . 
The resulting measure, named as Gaussian EoF in Ref. 11711 . 
is an upper bound to the true EoF and coincides with it for 
symmetric two-mode Gaussian states. 
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In general, we can define a Gaussian EM Ge as follows. 
For any pure Gaussian state ip with CM er p , one has 



G E {a p ) = E^), 



(22) 



where E can be any proper measure of entanglement of pure 
states, defined as a monotonically increasing function of the 
entropy of entanglement {i.e. the von Neumann entropy of the 
reduced density matrix of one party). 

For any mixed Gaussian state g with CM er, one has ifTvtl 



G e {(t)= inf G E (0- 

(T p <er 



(23) 



If the function E is taken to be exactly the entropy of entan- 
glement, then the corresponding Gaussian EM is known as 
Gaussian EoF fl7ll . In Ref. 13711 the properties of the Gaus- 
sian EoF have been further investigated, and interesting con- 
nections with the capacity of bosonic Gaussian channels have 
been established. 

In general, the definition Eq. d23i involves an optimization 
over all pure Gaussian states with CM er p smaller than the 
CM er of the mixed state whose entanglement one wishes 
to compute. Despite being a simpler optimization problem 
than that appearing in the definition Eq. ( 1191 of the true EoF 
(which, in CV systems, would imply considering decompo- 
sitions over all, Gaussian and non-Gaussian pure states), the 
Gaussian EMs cannot be expressed in a simple closed form, 
not even in the simplest instance of (nonsymmetric) two- 
mode Gaussian states. It is the aim of the present paper to 
compute Gaussian EMs for two relevant classes of, generally 
nonsymmetric, two-mode Gaussian states, namely the states 
of extremal (maximal and minimal) negativity at fixed global 
and local purities fl9l Eflll . which will be reviewed in Sec. 
IIV Al This will provide an insight into the problem of the or- 
dering 1 38] of two-mode Gaussian states with respect to dif- 
ferent measures of entanglement, leading to results somehow 
similar to those obtained for systems of two qubits 1 39], where 
in general the EoF and the negativity are found to be inequiv- 
alent. 

Before moving on to the explicit computations, let us re- 
call, as an important side remark, that any Gaussian EM is 
an entanglement monotone under Gaussian LOCC. The proof 
given in Sec. IV of Ref. 111711 for the Gaussian EoF, in fact, au- 
tomatically extends to every Gaussian EM constructed via the 
Gaussian convex roof of any proper measure E of pure-state 
entanglement. 



IV. GAUSSIAN ENTANGLEMENT MEASURES 
FOR TWO-MODE GAUSSIAN STATES 

The problem of evaluating Gaussian EMs for a generic 
two-mode Gaussian state has been solved in Ref. 0. How- 
ever, the explicit result contains so "cumbersome" expressions 
(involving the solutions of a fourth-order algebraic equation), 
that the authors of Ref. ifTvtl considered them not particularly 
useful to be reported explicitely in their paper. 

We recall here the computation procedure 11711 that we will 
need in the following. For any two-mode Gaussian state with 



CM er = <r s f in standard form Eq. (|9}, a generic Gaussian 
EM Ge is given by the entanglement E of the least entangled 
pure state with CM er p < er. Denoting by j q (respectively 
7 P ) the 2x2 submatrix obtained from er by canceling the 
even (resp. odd) rows and columns, we have, explicitely 



lq = 



a c + 
c+ b 



1p = 



a c_ 
c_ b 



(24) 



All the covariances relative to the "position" operators of the 
two modes are grouped in j q , and analogously for the "mo- 
mentum" operators in 7 p . The total CM can then be written as 
a direct sum er = j q © j p . Similarly, the CM of a generic pure 
two-mode Gaussian state in standard form (it has been proven 
that the CM of the optimal pure state has to be in standard 
form as well ifTvtl > can be written as er p = 7 P 7 P , where 
the global purity of the state imposes (7 P ) _1 = 7 P = T. The 
pure states involved in the definition of the Gaussian EM must 
thus fulfill the condition 



7p 1 < r < lq . 



(25) 



This problem is endowed with a nice geometric description 
JT7I1 . Writing the matrix T in the basis constituted by the iden- 
tity matrix and the three Pauli matrices, 



r = 



x + £3 xi 
Xi x - x 3 



(26) 



the expansion coefficients (xo, x%, X3) play the role of space- 
time coordinates in a three-dimensional Minkowski space. In 
this picture, for example, the rightmost inequality in Eq. (1251 
is satisfied by matrices T lying on a cone, which is equivalent 
to the (backwards) light cone of C q in the Minkowski space; 
and similarly for the leftmost inequality. Indeed, one can show 
that, for the optimal pure state cr p pt realizing the minimum in 
Eq. d23i . the two inequalities in Eq. d25i have to be simul- 
taneously saturated 1 17]. From a geometrical point of view, 
the optimal V has then to be found on the rim of the intersec- 
tion of the forward and the backward cones of 7~ 1 and j q , 
respectively. This is an ellipse, and one is left with the task 
of minimizing the entanglement E of er p = T (see 
Eq. (E3) for T lying on this ellipse lEfoll . 

At this point, let us pause to briefly recall that any pure 
two-mode Gaussian state er p is locally equivalent to a two- 
mode squeezed state with squeezing parameter r, described 
by a CM 



sq 



( cosh(2r) 


sinh(2r) 

V 





cosh(2r) 


— sinh(2r) 



sinh(2r) 


cosh(2r) 






— sinh(2r) 


cosh(2r) 



(27) 

The following statements are then equivalent: (i) E is a mono- 
tonically increasing function of the entropy of entanglement; 
(ii) E is a monotonically increasing function of the single- 
mode determinant m = Dct a. = Det (3 (see Eq. (|8}); (iii) 
E is a monotonically decreasing function of the local purity 
Hi = Hi = [i 2 (see Eq. (0); (iv) E is a monotonically de- 
creasing function of the smallest symplectic eigenvalue v p of 
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the partially transposed CM <x p ; (v) E is a monotonically in- 
creasing function of the squeezing parameter r. This chain 
of equivalences is immediately proven by simply recalling 
that a pure state is completely specified by its single-mode 
marginals, and that for a single-mode Gaussian state there 
is a unique symplectic invariant (the determinant), so that all 
conceivable entropic quantities are monotonically increasing 
functions of this invariant 1 20 ] . In particular, statement (ii) 
allows us to minimize directly the single-mode determinant 
over the ellipse: 



1 



Dctr 



(28) 



with r given by Eq. J26I . 

To simplify the calculations, one can move to the plane of 
the ellipse with a Lorentz boost which preserves the relations 
between all the cones; one can then choose the transformation 
so that the ellipse degenerates into a circle (with fixed radius), 
and introduce polar coordinates on this circle. The calcula- 



J 



tion of the Gaussian EM for any two-mode Gaussian state is 
thus finally reduced to the minimization of m from Eq. d28t . 
at given standard-form covariances of <x, as a function of the 
polar angle 9 on the circle 14011 . So far, this technique has 
been applied to the computation of the Gaussian EoF by min- 
imizing Eq. d28i numerically fPUl (see also El). In addition 
to that, as already mentioned, the Gaussian EoF has been ex- 
actly computed for symmetric states, and it has been proven 
that in this case the Gaussian EoF is the true EoF Q(J. 

In this work we present new analytical calculations of the 
Gaussian EMs for two relevant classes of nonsymmetric two- 
mode Gaussian states: the states of extremal negativities at 
fixed global and local purities 1 20], which will be introduced 
in the next subsection. We begin by writing the general ex- 
pression of the single-mode determinant Eq. d28l > in terms of 
the covariances of a generic two-mode state (see Eq. (|9jl) and 
of the polar angle 9. After some tedious but straightforward 
algebra, one finds 



mg(a,b, c+, c_) = 1 + < c+(ab — c 2 .) — c_ + cos 0d [a — b(ab — c?_)] [6 — a(ab — c 2 ,)] 



2 (ab ~ (?_) (a 2 + b 2 + 2c+c+) 

cos 9 [2abc 3 _ + (a 2 + b 2 ) c + c 2 _ + ((l - 26 2 ) a 2 + b 2 ) c_ - ab (a 2 + b 2 - 2) c+] 
^ [a - b(ab - c 2 _)} [b - a(ab - c 2 _)] 



sin 9 (a 2 - b 2 ) 



1 - 



[c+ (ab - c?_) +c- 



\ [a~ b(ab - c 2 _)} [b ~ a{ab - <?_)} 



(29) 



where we have assumed c+ > |c_| without any loss of gen- 
erality. This implies that, for any entangled state, c+ > and 
c_ < 0. The Gaussian EM (defined in terms of the function E 
on pure states, see Eq. J22l i) of a generic two-mode Gaussian 
state coincides then with the entanglement E computed on the 
pure state with m — m opt , with m op t = ming(me). Accord- 
ingly, the symplectic eigenvalue V- of the partial transpose of 
the corresponding optimal pure-state CM cr^ pt , realizing the 
infimum in Eq. j23l >. would read (see Eq. Eq. dl3l l) 

V-opt = v- (crf pt ) = ytm opt - y/m op t - 1 . (30) 

As an example, for the Gaussian EoF one has 

G Ef ((t) = h (v-aptimopt)) , (31) 

with h(x) defined by Eq. (12 1 1 . 

Finding the minimum of Eq. d29i analytically for a generic 
state is a difficult task. Numerical investigations show that 
the equation dgmg = can have from one to four physical 
solutions (in a period) corresponding to extremal points, and 



the global minimum can be attained in any of them depending 
on the parameters of the CM er under inspection. However, a 
closed solution can be found for two important classes of non- 
symmetric two-mode Gaussian states, as we will now show. 



A. Parametrization of two-mode covariance matrices and 
definition of extremal states 

We have shown in Refs. ifT^ Eoll that, at fixed global pu- 
rity /i = Tr q 2 of the global state g, and at fixed local purities 
fii t 2 = Tr g 2 2 of each of the two reduced single-mode states 
Qi = Tij^g, the smallest symplectic eigenvalue i>_ of the 
partial transpose of the CM a of a generic two-mode Gaus- 
sian state (which qualifies its separability by the PPT criterion, 
and quantifies its entanglement in terms of the negativities) is 
strictly bounded from above and from below. This entails the 
existence of two disjoint classes of extremal states, namely 
the states of maximum negativity for fixed global and local 
purities (GMEMS), and the states of minimum negativity for 
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fixed global and local purities (GLEMS) 1 20] . The negativi- 
ties of the two extremal classes of Gaussian states, moreover, 
turn out to remain very close to each other for all the possi- 
ble assignments of the three purities, allowing for a reliable 
experimental estimate of the negativity of a generic two-mode 
Gaussian state in terms of the average negativity lfl9ll . The 
latter is determined by knowledge of the three purities alone, 



which, in turn, may be experimentally measured in direct, pos- 
sibly efficient, ways lE^ll . 

Recalling these results, one can provide a very useful and 
insightful parametrization of the entangled two-mode Gaus- 
sian states in standard form (see also I43l0 . In fact, the co- 
efficients appearing in Eq. can be rewritten, in general, 
according to the following, useful parametrization: 



a = s + d , b = s — d , 



(32) 



c± 



4Vs 2 - d? 



± 



-1 2 



Ad? + -{g* + l){\-l)-{2d? + g) (A + l) 



4s 2 + _ ( 5 2 + 1) (A - 1) - {2d 2 + g) (A + 1) 



(33) 



r 



where the two local purities are regulated by the parameters 
s and d, being fii = (s + fi2 = [s — <2) _1 , and the 

global purity is n = g~ x . The coefficient A embodies the only 
remaining degree of freedom needed for the complete deter- 
mination of the negativities, once the three purities have been 
fixed. It ranges from the minimum A = — 1 (corresponding to 
the GLEMS) to the maximum A = +1 (corresponding to the 
GMEMS). Therefore, as it varies, A encompasses all possible 
entangled two-mode Gaussian states compatible with a given 
set of assigned values of the purities. The constraints that the 
parameters s, d, g must obey for Eq. © to denote a proper 
CM of a physical state are: s > 1, \d\ < s — 1, and 



S >2|d| + 1, 



(34) 



If the global purity is large enough so that Ineq. ( 1341 is 
saturated, GMEMS and GLEMS coincide, the CM becomes 
independent of A, and the two classes of extremal states co- 
alesce into a unique class, completely determined by the 
marginals s and d. We denote these states as GMEMMS 
l2(ill . that is, Gaussian two-mode states of maximal negativity 
at fixed local purities. Their CM is simply characterized by 
c± = ±^/s 2 — (d + l) 2 , where we have assumed, without 
any loss of generality, that d > (corresponding to choose, 
for instance, mode 1 as the more mixed one: \i\ < fi2). 

In general 1 19], a GMEMS (A = +1) is entangled for 



9 < 2s 



(35) 



while a GLEMS (A = —1) is entangled for smaller g, namely 

g< y/2(s* + d?) - 1 . (36) 

To have a physical insight on these peculiar two-mode states, 
let us recall 1 20] that GMEMS are simply nonsymmetric ther- 
mal squeezed states, usually referred to as maximally entan- 
gled mixed states in CV systems. On the other hand, GLEMS 
are mixed states of partial minimum uncertainty, in the sense 



that the smallest symplectic eigenvalue of their CM is equal 
to 1, saturating the uncertainty inequality (II 11 1. 

We are now equipped with the necessary tools, and in the 
next subsection we move on to compute Gaussian EMs for 
the two extremal classes of nonsymmetric two-mode Gaus- 
sian states, the GLEMS and the GMEMS. 



Gaussian entanglement of minimum-negativity states 
(GLEMS) 



We want to find the optimal pure state cr^ pt entering in the 
definition Eq. d23l of the Gaussian EM. To do this, we have 
to minimize the single-mode determinant of cr^ pt , given by 
Eq. (I29> . over the angle 9. It turns out that, for a generic 
GLEMS, the coefficient of sin# in the last line of Eq. d29t 
vanishes, and the expression of the single-mode determinant 
reduces to the simplified form 



.CLEMS 



[AcosO + Bf 



2(ab-c 2 _)[{g 2 -\)cos6 + g 2 + l} ' 



(37) 



with A — c + (ab — c 2 _) + c_, B = c + (ab — c 2 _) — c_, 
and a, b, c± are the covariances of GLEMS, obtained from 
Eqs. I32l33i setting A = -1. 

The only relevant solutions (excluding the unphysical and 
the trivial ones) of the equation demo = are 9 = ir and 



9 = ±9* 



3 + 5 2 



2c_ 



1 



■r 



c+(ab — c_) + c- 



Studying the second derivative d 2 mo for 9 
immediately that, for 



> 4 - 



2c+(ab-c 2 



ir one finds 



(38) 
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(remember that c_ < 0), the solution 9 = n is a minimum. In 
this range of parameters, the other solution 8 = 9* is unphys- 
ical (in fact | cos 9* \ > 1), so me =7r is the global minimum. 
When, instead, Ineq. J38l > is violated, mg has a local maxi- 
mum for 9 = 7r and two minima appear at 9 = ±9*. The 



global minimum is attained in any of the two, given that, for 
GLEMS, rag is invariant under reflection with respect to the 
axis 9 = it. Collecting, substituting, and simplifying the ob- 
tained expressions, we arrive at the final result for the optimal 

/7t: 



J 



''opt 



16s 2 rf 2 



-g 4 +2(2d 2 +2s 2 + l)g 2 -(4d 2 -l)(4s 2 -l)-\/5 
s7^ ' 



g > y/2(s 2 + <P) — 1 [separable state] 



^ (4^ + i) d2+ , 2 +4,ypTw^| d | £ g < V 2(g2 + d2) 



2\d\ + l<g< 



(4s 2 + l)d 2 +s 2 +4s\d\^{s 2 + l)d 2 +s 2 
d 2 +s 2 



_£39) 



Here 8 = (2d - g - l)(2d -g + I) (2d + g- l)(2d + g + 
l)(g - 2s - l)(g - 2s + l)(g + 2s - l)(g + 2s + 1). 

Immediate inspection crucially reveals that m'^ ms is not 
in general a function of the symplectic eigenvalue i>_ alone. 
Therefore, unfortunately, the Gaussian EMs, and in particular, 
the Gaussian EoF, are not equivalent to the negativities for 
GLEMS. Further remarks will be given in the following, when 
the Gaussian EMs of GLEMS and GMEMS will be compared 
and their relationship with the negativities will be elucidated. 



p± _ p± 



turns out 



C. Gaussian entanglement of maximum-negativity states 
(GMEMS) 

The minimization of nig from Eq. ( I29> can be carried out 
in a simpler way in the case of GMEMS, whose covari- 
ances can be retrieved from Eq. J33l > setting A = 1. First 
of all, one can notice that, when expressed as a function of 
the Minkowski coordinates (xq, xx, £3), corresponding to the 
submatrix T Eq. d26t of the pure state er p = r®L ] en- 
tering in the optimization problem Eq. (1231 . the single-mode 
determimant m of <r p is globally minimized for £3 = 0. In 
fact, from Eq. ( 1281 . m is minimal, with respect to X3, when 
Det F = — x\ — x\ is maximal. Next, one can show that 
for GMEMS there always exists a matrix T, with X3 = 0, 
which is a simultaneous solution of the two matrix equations 
obtained by imposing the saturation of the two sides of in- 
equality ( I25> . As a consequence of the above discussion, this 
matrix would denote the optimal pure state cr p pt . Solving the 
system of equations Det (7^ — T) = Det (T — jp 1 ) = 0, 
where the matrices involved are explicitely defined combin- 
ing Eq. (I24> and Eq. (1331 with A = 1, one finds the following 
two solutions for the coordinates xo and x\ : 



± _ (9 + i)g ± Vl(g - i) 2 - 4rf 2 ] M 2 + g - .9) 

C ° 2(d?+g) 
± (g + l)V-d 2 + s 2 - g ± Sy /(g - l) 2 - 4d 2 
Cl 2 (d 2 + g) 



(40) 



The corresponding pure state a 
to be, in both cases, a two-mode squeezed state described by 
a CM of the form Eq. 1271 . with cosh(2r) = x J . Because the 
single-mode determinant 777 = cosh 2 (2r) for these states, the 
optimal 777 for GMEMS is simply equal to (a;^) 2 . Summariz- 
ing, 



1, 



g > 2s — 1 [separable state] ; 



-.GMEMS 

l opt 



{(g+l)s-^/l(g-l) 2 -4d 2 ](-d 2 + s 2 -g)Y 
4(d 2 +g) 2 

2\d\ + 1 < g < 2s - 1 . 



(41) 

Once again, also for the class of GMEMS the Gaussian EMs 
are not simple functions of the symplectic eigenvalue V- 
alone. Consequently, they provide a quantification of CV en- 
tanglement of GMEMS inequivalent to the one determined by 
the negativities. Furthermore, we will now show how these re- 
sult raise the problem of the ordering of two-mode Gaussian 
states according to their degree of entanglement, as quantified 
by different families of entanglement measures. 



V. EXTREMAL ORDERING OF 

TWO-MODE GAUSSIAN STATES 

Entanglement is a physical quantity. It has a definite math- 
ematical origin within the framework of quantum mechanics, 
and its conceptual meaning in the end stems from and is rooted 
in the existence of the superposition principle. Further, en- 
tanglement has a fundamental operative interpretation as that 
resource that in principle enables information processing and 
communication in better-than-classical realizations |8]. One 
would then expect that, picking two states qa and qb out of 
a certain (subset of) Hilbert space, the question "Is qa more 
entangled than gs?" should have a unique, well-defined an- 
swer, independent of the measure that one chooses to quantify 
entanglement. But, contrary to the common expectations, this 
is generally not the case for mixed states. Different measures 




FIG. 1: (color online). Comparison between the ordering induced by 
Gaussian EMs on the classes of states with extremal (maximal and 
minimal) negativities. This extremal ordering of the set of entangled 
two-mode Gaussian states is studied in the space of the CM's param- 
eters {s, d, g}, related to the global and local purities by the relations 
[Ax = (s + H2 = (s — d)" 1 and fi = g~ x . The intermediate, 

meshed surface is constituted by those global and local mixednesses 
such that the Gaussian EMs give equal values for the correspond- 
ing GMEMS (states of maximal negativities) and GLEMS (states of 
minimal negativities). Below this surface, the extremal ordering is 
inverted (GMEMS have less Gaussian EM than GLEMS). Above it, 
the extremal ordering is preserved (GMEMS have more Gaussian 
EM than GLEMS). However, it must be noted that this does not ex- 
clude that the individual orderings induced by the negativities and by 
the Gaussian EMs on a pair of non-extremal states may still be in- 
verted in this region. Above the uppermost, lighter surface, GLEMS 
are separable states, so that the extremal ordering is trivially pre- 
served. Below the lowermost, darker surface, no physical two-mode 
Gaussian states can exist. All the quantities plotted are dimension- 
less. 



of entanglement will in general induce different, inequivalent 
orderings on the set of entangled states belonging to a given 
Hilbert space 13811 . as they usually measure different aspects 
of quantum correlations existing in generic mixed states. 

In the context of C V systems, when one restricts to symmet- 
ric, two-mode Gaussian states, which include all pure states, 
the known computable measures of entanglement all correctly 
induce the same ordering on the set of entangled states. We 
will now show that, indeed, this nice feature is not preserved 
moving to mixed, nonsymmetric two-mode Gaussian states. 
We aim at comparing Gaussian EMs and negativities on the 
two extremal classes of two-mode Gaussian states 1 20], intro- 
ducing thus the concept of extremal ordering. At fixed global 
and local purities, the negativity of GMEMS (which is the 
maximal one) is obviously always greater than the negativity 
of GLEMS (which is the minimal one). If for the same val- 
ues of purities the Gaussian EMs of GMEMS are larger than 



II " 




preserved 
extremal 
ordering 



1 2 3 4 5 

b = fi 2 - 1 

FIG. 2: (color online). Summary of entanglement properties of two- 
mode Gaussian states, in the projected space of the local mixedness 
b — fi 2 1 of mode 2, and of the global mixedness g — [lT 1 , while 
the local mixedness of mode 1 is kept fixed at a reference value 
a — fj,^ 1 — 5. Below the thick curve, obtained imposing the equal- 
ity in Ineq. 1421 . the Gaussian EMs yield GLEMS more entangled 
than GMEMS, at fixed purities: the extremal ordering is thus in- 
verted. Above the thick curve, the extremal ordering is preserved. In 
the coexistence region (see Ref. tl^l ). GMEMS are entangled while 
GLEMS are separable. The boundaries of this region are given by 
Eq. I36i (dashed line) and Eq. J35I (dash-dotted line). In the separa- 
bility region, GMEMS are separable too, so all two-mode Gaussian 
states whose purities lie in that region are not entangled. The shaded 
regions cannot contain any phisical two-mode Gaussian state. All the 
quantities plotted are dimensionless. 



those of GLEMS, we will say that the extremal ordering is 
preserved. Otherwise, the extremal ordering is inverted. In 
this latter case, which is clearly the most intriguing, the states 
of minimal negativities are more entangled, with respect to 
Gaussian EMs, than the states of maximal negativities, and 
the inequivalence of the orderings, induced by the two differ- 
ent families of entanglement measures, becomes manifest. 



The problem can be easily stated. By comparing m 



from Eq. d39l and 



GLEMS 

opt 



i opt 



from Eq. J41L one has that in the 



range of global and local purities, or, equivalently, of parame- 
ters {s, d, <?}, such that 



"opt 



> m 



opt 



(42) 



the extremal ordering is preserved. When Ineq. i42l is vi- 
olated, the extremal ordering is inverted. The boundary be- 
tween the two regions, which can be found imposing the 



equality m 



GMEMS 

opt 



l opt 



\ yields the range of global 



and local purities such that the corresponding GMEMS and 
GLEMS, despite having different negativities, have equal 
Gaussian EMs. This boundary surface can be found numer- 
ically, and the result is shown in the 3D plot of Fig. ^ 

One can see, as a crucial result, that a region where the 
extremal ordering is inverted does indeed exist. The Gaus- 
sian EMs and the negativities are thus definitely not equiva- 
lent for the quantification of entanglement in nonsymmetric 
two-mode Gaussian states. The interpretation of this result 
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is quite puzzling. On the one hand, one could think that the 
ordering induced by the negativities is a natural one, due to 
the fact that such measures of entanglement are directly in- 
spired by the necessary and sufficient PPT criterion for sep- 
arability. Thus, one would expect that the ordering induced 
by the negativities should be preserved by any bona fide mea- 
sure of entanglement, especially if one considers that the ex- 
tremal states, GLEMS and GMEMS, have a clear physical in- 
terpretation [20]. Therefore, as the Gaussian EoF is an upper 
bound to the true EoF, one could be tempted to take this result 
as an evidence that the Gaussian EoF overestimates the true 
EoF, at least for GLEMS, and that, moreover, the true EoF 
of GLEMS should be lower than the true EoF of GMEMS, 
at fixed values of the purities. If this were the case, the true 
EoF would not coincide with the Gaussian EoF, whose eval- 
uation would consequently necessarily involve a decomposi- 
tion over non-Gaussian states. However, this is only a qual- 
itative/speculative argument: proving or disproving that the 
Gaussian EoF is the true EoF for any two-mode Gaussian 
state is still an open question under lively debate lll8ll . 

On the other hand, one could take the simplest discrete- 
variable instance, constituted by a two-qubit system, as a 
test-case for comparison. There, although for pure states the 
negativity coincides with the concurrence, an entanglement 
monotone equivalent to the EoF for all states of two qubits 
J44ll . the two measures cease to be equivalent for mixed states, 
and the orderings they induce on the set of entangled states 
can be different |39]. This analogy seems to support again 
the stand that, in the arena of mixed states, a unique mea- 
sure of entanglement is a chimera and cannot really be ex- 
pected, due to the different operative meanings and physical 
processes (in the cases when it has been possible to identify 
them) that are associated to each definition: one could think, 
for instance, of the operative difference existing between the 
definitions of distillable entanglement and entanglement cost. 
In other words, from this point of view, each inequivalent 
measure of entanglement introduced for mixed states should 
capture physically distinct aspects of quantum correlations ex- 
isting in these states. Then, joining this kind of outlook, one 
could hope that the Gaussian EMs might still be considered as 
proper measures of CV entanglement, especially if one were 
able to prove the conjecture that the Gaussian EoF is the true 
EoF for a broader class of Gaussian states beyond the symmet- 
ric ones. One could then live on with the existence of inverted 
orderings of entangled states, and see it as a not so annoying 
problem. 

Whatever be the case, we have shown that two different 
families of measures of CV entanglement can induce differ- 
ent orderings on the set of two-mode entangled states. This 
is more clearly illustrated in Fig. |3] where we keep fixed one 
of the local mixednesses and we classify, in the space of the 
other local mixedness and of the global mixedness, the differ- 
ent regions related to entanglement and extremal ordering of 
two-mode Gaussian states, improving and completing a sim- 
ilar diagram previously introduced in Ref. 1 19] to describe 
separability in the space of purities. 



VI. GAUSSIAN MEASURES OF ENTANGLEMENT 
VERSUS NEGATIVITIES 

In this section we wish to give a more direct comparison 
of the two families of entanglement measures for two-mode 
Gaussian states. In particular, we are interested in finding the 
maximum and minimum values of one of the two measures, if 
the other is kept fixed. A very similar analysis has been per- 
formed by Verstraete et al. 1 39], in their comparative analysis 
of the negativity and the concurrence for states of two-qubit 
systems. 

Here it is useful to perform the comparison directly between 
the symplectic eigenvalue V-{cr) of the partially transposed 
CM cr of a generic two-mode Gaussian state with CM cr, 
and the symplectic eigenvalue i>_ (cr£, 4 ) of the partially trans- 
posed CM cr^ t of the optimal pure state with CM cr^ pt , which 



minimizes Eq. (I23> . In fact, the negativities are all monotoni- 
cally decreasing functions of (cr), while the Gaussian EMs 
are all monotonically decreasing functions of j>_(er^ t ). 

To start with, let us recall once more that for pure states and 
for mixed symmetric states (in the set of two-mode Gaussian 
states), the two quantities coincide. For nonsymmetric states, 
one can immediately prove the following bound 



(a- opt) < v-{(r) . 



(43) 



In fact, from Eq. (1231 . cr^L 



< cr fl7ll . For positive matrices, 
A > B implies > b^, where the a^-s (resp. bkS) denote the 
ordered symplectic eigenvalues of A (resp. B) B45I1 . Because 
the ordering A > B is preserved under partial transposition, 
Ineq. (I43i holds true. This fact induces a characterization of 
symmetric states, which saturate Ineq. ( I43> . as the two-mode 
Gaussian states with minimal Gaussian EMs at fixed negativ- 
ities. 

It is then natural to raise the question whether an upper 
bound on the Gaussian EMs at fixed negativities exists as well. 
It seems hard to address this question directly, as one lacks 
a closed expression for the Gaussian EMs of generic states. 
But we can promptly give partial answers if we restrict to the 
classes of GLEMS and of GMEMS, for which the Gaussian 
EMs have been explicitely computed in the previous section. 

Let us begin with the GLEMS. We can compute 
the squared symplectic eigenvalue i} 2 (er GLEMS ) = 



4(s 2 + d 2 ) - g 2 - 1 - y(4(.s 2 + d 2 ) - g 2 - \) A - 4g 2 

Next, we can reparametrize the CM (obtained by Eq. J33I 
with A = —1) to make i>_ appear explicitely, namely 

yJvl[A(3 2 +Cp) 



/2- 
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1 



'-]/(! 



[ ). At this point, 



one can study the piecewise function m^ f EMS from Eq. ( I39i . 
and find out that it is a convex function of d in the whole 
space of parameters corresponding to entangled states. 
Hence, TOgp t BMS , and thus the Gaussian EM, is maximized at 
the boundary \d\ = (2us — i^L — l)/2, resulting from the 
saturation of Ineq. ( I34l >. The states maximizing Gaussian EMs 
at fixed negativities, if we restrict to the class of GLEMS, 
have then to be found in the subclass of GMEMMS (states 
of maximal negativity for fixed marginals 1 20], defined after 
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FIG. 3: (color online). Comparison between Gaussian EMs and neg- 
ativities for two-mode Gaussian states. On the horizontal axis we 
plot the symplectic eigenvalue v~ (cr) of the partially transposed CM 
ct of a generic two-mode Gaussian state with CM cr. On the verti- 
cal axis we plot the symplectic eigenvalue i/_(cr^, t ) of the partially 
transposed CM cr^ p t of the optimal pure state with CM cr^ pt , which 
minimizes Eq. I23t . The negativities are all monotonically decreas- 
ing functions of j/_(cr), while the Gaussian EMs are all monoton- 
ically decreasing functions of i>-(cr^ pt ). The equation of the two 
boundary curves are obtained from the saturation of Ineq. 1431 (upper 
bound) and Ineq. 1461 (lower bound), respectively. The dots repre- 
sent 50 000 randomly generated CMs of two-mode Gaussian states. 
Of up to 1 million random CMs, none has been found to lie below 
the lower solid-line curve, enforcing the conjecture that it be an ab- 
solute boundary for all two-mode Gaussian states. All the quantities 
plotted are dimensionless. 

Ineq. 034H. depending on the parameter s and on the eigen- 
value itself, which completely determines the negativity). 
For these states, 

™? P T MMS (-,*-) = ( i3^T2^ ) ■ (44) 

The further optimization over s is straightforward because 
TO opt EMMS i s an increasing function of s, so its global max- 
imum is attained for s — > oo. In this limit, one has simply 

m GMEMM S (£>_) = _L. (45) 

From Eq. (|30j, one thus finds that for all GLEMS the follow- 
ing bound holds 

v- « t ) > (l - . (46) 

One can of course perform a similar analysis for GMEMS. 
But, after analogous reasonings and computations, what one 
finds is exactly the same result. This is not so surprising, keep- 
ing in mind that GMEMS, GLEMS and all two-mode Gaus- 
sian states with generic s and d but with global mixedness 
g saturating Ineq. j34t . collapse into the same family of two- 
mode Gaussian states, the GMEMMS, completely determined 
by the local single-mode properties (they can be viewed as a 
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FIG. 4: (color online). Comparison between the Gaussian entangle- 
ment of formation Ge f and the logarithmic negativity £V for two- 
mode Gaussian states. Symmetric states accomodate on the lower 
boundary (solid line), determined by the saturation of Ineq. 1471 . 
GMEMMS with infinite, average local mixedness, lie on the dashed 
line, whose defining equation is obtained from the saturation of 
Ineq. |48}. All GMEMS and GLEMS lie below the dashed line. The 
latter is conjectured, with strong numerical support, to be the upper 
boundary for the Gaussian EoF of all two-mode Gaussian states, at 
fixed negativity. All the quantities plotted are dimensionless. 



generalization of the pure two-mode states: the symmetric 
GMEMMS are in fact pure). Hence, the bound of Ineq. J46l . 
limiting the Gaussian EMs from above at fixed negativities, 
must hold for all GMEMS as well. 

At this point, it is tempting to conjecture that Ineq. J46l > 
holds for all two-mode Gaussian states. Unfortunately, the 
lack of a closed, simple expression for the Gaussian EM of 
a generic state makes the proof of this conjecture impossible, 
at the present time. However, one can show, by analytical 
power-series expansions of Eq. i29\ . truncated to the leading 
order in the infinitesimal increments, that, for any infinitesi- 
mal variation of the parameters of a generic CM around the 
limiting values characterizing GMEMMS, the Gaussian EMs 
of the resulting states lie always below the boundary imposed 
by the corresponding GMEMMS with the same v_. In this 
sense, the GMEMMS are, at least, a local maximum for the 
Gaussian EM versus negativity problem. Furthermore, ex- 
tensive numerical investigations of up to a million CMs of 
randomly generated two-mode Gaussian states, provide con- 
firmatory evidence that GMEMMS attain indeed the global 
maximum (see Fig. [3}. We can thus quite confidently con- 
jecture, however, at the moment, without a complete formal 
proof of the statement, that GMEMMS, in the limit of infinite 
average local mixedness (s — > oo), are the states of maxi- 
mal Gaussian EMs at fixed negativities, among all two-mode 
Gaussian states. 

A direct comparison between the two prototypical represen- 
tatives of the two families of entanglement measures, respec- 
tively the Gaussian EoF Ge f and the logarithmic negativity 
Etf, is plotted in Fig. |4] For any fixed value of E^f, Ineq. i43i 
provides in fact a rigorous lower bound on Ge f , namely 

G Ef > h[exp(-E„)] , (47) 
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while Ineq. J46> provides the conjectured lower bound 

G Ef < h \exp{E M ) (l - v/1 - exp(-2S AA ) 



(48) 



where we exploited Eqs. ( 11813 11 and h[x] is given by Eq. i !2 It . 

The existence of lower and upper bounds on the Gaussian 
EMs at fixed negativities (the latter strictly proven only for 
extremal states), limits to some extent the inequivalence aris- 
ing between the two families of entanglement measures, for 
nonsymmetric two-mode Gaussian states. 



VII. SUMMARY AND OUTLOOK 

In this work we focused on the simplest conceivable states 
of a bipartite CV system: two-mode Gaussian states. We have 
shown that, even in this simple instance, the theory of quan- 
tum entanglement hides several subtleties and reveals some 
surprising aspects. In particular, we have studied the rela- 
tions existing between different computable measures of en- 
tanglement, showing how the negativities (including the stan- 
dard logarithmic negativity) and the Gaussian convex-roof ex- 
tended measures (Gaussian EMs, including the Gaussian en- 
tanglement of formation 1 17]) are inequivalent entanglement 
quantificators for nonsymmetric two-mode Gaussian states. 
We have computed Gaussian EMs explicitely for the two 
classes of two-mode Gaussian states having extremal (max- 
imal and minimal) negativities at fixed purities 1 20] . We have 
highlighted how, in a certain range of values of the global and 
local purities, the ordering on the set of entangled states, as 
induced by the Gaussian EMs, is inverted with respect to that 
induced by the negativities. The question whether a certain 
Gaussian state is more entangled than another, thus, has no 
definite answer, not even when only extremal states are con- 
sidered, as the answer comes to depend on the measure of en- 
tanglement one chooses. Extended comments on the possible 
meanings and consequences of the existence of inequivalente 
orderings of entangled states have been given in SectionMand 
in Section IVT1 Furthermore, we have proven the existence of 
a lower bound holding for the Gaussian EMs at fixed negativ- 
ities, and that this bound is saturated by two-mode symmet- 
ric Gaussian states. Finally, we have provided some strong 
numerical evidence, and partial analytical proofs restricted to 



extremal states, that an upper bound on the Gaussian EMs at 
fixed negativities exists as well, and is saturated by states of 
maximal negativity for given marginals, in the limit of infinite 
average local mixedness. 

We believe that our results will raise renewed interest in the 
problem of the quantification of entanglement in CV systems, 
which seemed fairly well understood in the special instance 
of two-mode Gaussian states. Moreover, we hope that the 
present work may constitute a first step toward the solution of 
more general problems concerning the entanglement of Gaus- 
sian states, such as the computation of the entanglement of 
formation for generic two-mode Gaussian states II 811 . and the 
proof of its identity with the Gaussian EoF in a larger class of 
Gaussian states beyond the symmetric instance. On the other 
hand, the explicit expressions, computed in the present work, 
now available for the Gaussian EoF of GMEMS and GLEMS, 
might serve as well as a basis to find an explicit counterex- 
ample to the conjecture that the decomposition over all pure 
Gaussian states, in the definition of the EoF, is the optimal one 
for all two-mode Gaussian states. 

Finally, the results collected in the present work might 
prove useful as well in the task of quantifying multipartite 
entanglement of Gaussian states. For instance, we should 
mention here that any two-mode reduction of a pure three- 
mode Gaussian state is a GLEMS, as a consequence of the 
Schmidt decomposition operated at the CM level 1 46] . There- 
fore, thanks to the results that we have derived here, its Gaus- 
sian EoF can be explicitely computed, and can be compared 
with the entropy of entanglement between one reference mode 
and the remaining two in the global state. One has then avail- 
able the tools and can apply them to investigate the sharing 
structure of multipartite CV entanglement of three-mode, and, 
more generally, multimode Gaussian states l47ll . 
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